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Abstract 

We develop a rigorous treatment of discontinuous stochastic unitary 
evolution for a system of quantum particles that interacts singularly with 
quantum “bubbles” at random instants of time. This model of a ’’cloud 
chamber” allows to watch and follow with a quantum particle trajectory 
like in cloud chamber by sequential unsharp localization of spontaneous 
scatterings of the bubbles. Thus, the continuous reduction and sponta¬ 
neous localization theory is obtained as the result of quantum filtering 
theory, i.e., a theory describing the conditioning of the a priori quantum 
state by the measurement data. We show that in the case of indistinguish¬ 
able particles the a posteriori dynamics is mixing, giving rise to an irre¬ 
versible Boltzmann-type reduction equation. The latter coincides with the 
nonstochastic Schrodinger equation only in the mean field approximation, 
whereas the central limit yields Gaussian mixing fluctuations described 
by stochastic reduction equations of diffusive type. 


1 Introduction 

The quantum measurement theory based on the ordinary von Neumann reduc¬ 
tion postulate applies neither to instantaneous observations with continuous 
spectra nor to continual (continuous in time) measurements. Although such 
phenomena can be described in the more general framework of Ludwig’s Davies- 
Lewies operational approach n-ia, there is a particular interest in describing 
quantum measurements by concrete Hamiltonian models from which the oper¬ 
ational description can be derived by an averaging procedure. Perhaps the first 
model of such kind for instantaneous unsharp measurement of particle local¬ 
ization was given by von Neumann |^. He considered the singular interaction 

'^On leave of absence from Mathematics Department of Nottingham University, NG7 2RD, 
UK. 
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Hamiltonian 


K{t) = = ^ (1.1) 

1 dj/ 1^0, 

with Dirac (5-function, producing the translation operator 

h,(t)dt| (1.2) 

at time t = 0 of the measurement. Here x is the position of the particle and y is 
the pointer position given by the (^-coordinate of a quantum meter. The particle 
scattering operator S = {sa,} applied to the (generalized) position eigen-vectors 
lx) as S'|x) = |x)sx does not affect the position x of the particle but changes 
the meter coordinate q to y = s\.qsx = x -\- q. This implies that the initial wave 
function ip^{x,y) of the system “particle plus meter” is transformed into 

■ipix, y) = Sxipoix, y) = ip^ix, y - x). (1.3) 

If in the initial state the particle and the meter were not coherent, v) — 
q{x)fo{y), and the wave function /q of the meter was fixed, then one can obtain 
the unitary transformation S : tpQ s^ipQ via a family {F(y)} of reduction 
transformations F{y) : ij i—> f,{y)r] for the particle vector-states q. Specifically, 
Eq. (11.311 can be defined as 

y) = fo{y - x)q{x) = fx{y)q{x). (1.4) 

The linear nonunitary operators F[y) = s,fo{y) act on |x) as the multiplication 
Fiy)\x) = \x)fx{y) by fx{y) = fo{y - x) = Sxfo{y) and would give a sharp 
localization of any particle wave function q{x) at the point x = y ot the pointer 
position provided that the wave function foiy) could be initially localized at 
y = 0. But there are no such sharply localized quantum states for the continuous 
pointer, and the best that one can do is to take a wave packet foiy), say, of the 
Gaussian form 

/o(y) = expj - (1.5) 

This results in the unsharp localization 

= expj-^ (x - y)^|??(x) (1.6) 

of any particle wave function qix) about the observed value y of the pointer, 
normalized to the probability density 

Voiv) = j\hiy - x)\^\q{x)\^ X . (1.7) 

For commuting operators x and y, this is equivalent to the classical measurement 
model y = x + q for unsharp measurement of an unknown signal x via the sharp 
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measurement of the signal plus Gaussian noise q with given probability density 

l/o(<z)p. 

These arguments illustrate how to interpret the reduction model involving 
the continuous spectrum of a quantum measurement as a Hamiltonian inter¬ 
action model with nondemolition observation for a quantum object using the 
measurement of the pointer coordinate of the quantum meter. Since von Neu¬ 
mann introduced this approach, it was used by numerous other authors izmi 
for the derivation of a generalized reduction rj i—> F{y)rj that would replace the 
von Neumann postulate 77 1 —*■ E(jj)r] given by orthoprojections {E(y)}. As ex¬ 
tended to nondemolition observations continual in time this approach 

consists in using the quantum filtering method for the derivation of nonunitary 
stochastic wave equations describing the quantum dynamics under the observa¬ 
tion. Since a particular type of such equations has been taken as a postulate in 
the phenomenological theory of continuous reduction and spontaneous localiza¬ 
tion m-m, the question arises whether it is possible to obtain this equation 
from an appropriate Schrodinger equation. Here we shall show how this can be 
done by second quantization of the interaction Hamiltonian considered by von 
Neumann, obtaining a stochastic model of continual nondemolition observation 
for the position of a quantum particle by counting some other quanta. But 
first we show that even the projection postulate can be derived in the frame¬ 
work of this approach with the suggested Hamiltonian interaction and a proper 
nondemolition observation. 


2 Hamiltonian Reduction Model 


Let be a Hilbert space, called the state space of a particle, and let R = {R°‘ : 
a = l,...,c?} be selfadjoint operators in Ti with either integer or continuous 
spectrum or d G {1,2,3}. Let k > 0 be a scaling parameter. If the 
operators commute, one can regard the scaled vector operator kR as the position 
X of the particle in R'^ or in the d-dimensional lattice if it is quantized. In 
the x-representation, the operator acts as the multiplication kR°‘\x) = |x)A(a;“) 
by A(a:“) = k[x°‘/k\, where [xj G Z'^ denotes the integer part of the vector 
X = {a;“ : a = 1,..., d}. 

A quantum meter with continuous (A = R'^) or lattice (A = eZ"^) pointer 
scale is described by the Hilbert space L^(A) of complex-valued functions / : 
A ^ C square integrable in the sense that ||/|P = / |/(y)pdA < 00 , where 


Jf{y) dA = 


EyGA/(y)e'*. 

/yeA /(y)dy, 


if A = eZ'^, 
if A = R'^, 


dA = e^ 
dA = dy. 


Consider a moving particle with Hamiltonian H in Ji. Its singular evolution 
corresponding to the position measurement at time t = 0 is described in the 
product space 0.1=0.® L^(A) by the time-dependent Hamiltonian 


Hx{t) = Ho + Kli®5{t)V. 


( 2 . 1 ) 
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Here Hq = H ® 1, where 1 is the identity operator in L?‘{K), and R 0 P = 
'Y^PaR°‘, where Pa, a = are operator components of the meter mo¬ 

mentum vector P = (Pi,..., Pd) with the pointer coordinate vector q = {g“ : 
a = 1,... ,c?} given by the multiplication operators (q’“/)(y) = 2 /“/(y), y G A 
in the Hilbert space A^(A) of its quantum states such that 

d 

ly“/(y)PdA < oo. 

a—1 

The operators Pa = —ihd/dy^ have the bounded spectrum {'Kh)[—\/£, for 
A = eU^ and are defined by the matrix elements 

{y\Pc. |q) = JPaP^y-^^^/^dp , q G A , 

where (y|q) = ^(y — q), y G A, are the (generalized) eigenfunctions of the 
coordinate operators q°‘ normalized with respect to the measure dA = in the 
discrete case y°‘ G eZ. This generates the shift operators US in T^(A) with 
X G A and the scattering operator 

S, = exp(-rKR®P,)=|)- (2.2) 

in Til, where P* = ItP, It = 1 if t > 0, and It = 0 otherwise. 

The singular time dependence of the Hamiltonian eu makes it impossible 
to define the Schrodinger equation ihdijj/dt = in the usual sense. But 

one can define a unitary stochastic evolution Ui{t) : Hi —> Hi for some tg < 0 
as the single-jump unitary process 

Pi(t) =exp(^-^^ Hi{s)ds'^ = 

provided that [P“, H] = 0. Usually the commutator is not zero, and the unitary 
evolution corresponding to eu must be redefined in terms of the solution 
ip(t) = Ui{t)'ipQ to the regularized wave equation. This can be done in terms of 
the forward differentials 

d-ipd) = ili{t + dt) - dlt = It+dt - It; 

namely, we define the generalized Schrodinger equation as 

dipd) + ^HQil){t)dt ={S - I)'ijj{t)dlt, tpito) = ipg . (2.3) 

Proposition 1 For every ipQ G H the difference equation lH.f-Al has a unique 
solution corresponding to the initial value to ^ 0; this solution is given by the 
unitary operator 

Ui{t) = C/o(t-to)*S't(-to), where Uo{t) = eyip{-iHot/h}, St{r) = C/g(^)'5'tUo(r). 
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Proof. Let us rewrite the generalized Schrodinger equation in the integral 
form ^ 

^(t) = - I)i:{r)dlr'^ (2.4) 

(the equivalence of (12.411 to Il2.dll can be shown by straightforward differentia¬ 
tion). We can write ^/>(0) = {7o(—to)V^o ^he solution (12.411 when t = 0 since 
Jtg ^(^)dlr = 0 for any to ^ 0; hence tj}{t) can be rewritten as 

m = C/o(t)(^(0) + {S - l)lti^{0)) = Uoit)St^{0) 


for any to < 0 since 



1 *^( 0 ) 


and {S — I)lt = St — I by the definition of the scattering operator 12.211 . Thus, 
U(t) is the unitary operator 


Uoit)StUo{—to) — Uo{t — to)St{—to)- 

This gives the solution to equation (lOl for to = 0 as well, since tp{t) = 
Uo{t)StipiO) is equal to ^po t = 0. ■ 

The rescaled pointer qt = l^q switched on at the instant t = 0 of the 
scattering is described in the space Ti 0 L^(A) by the operators Q* = qt. 

In the Heisenberg picture Ul{t)QtUi{t) these operators are described for any 
to < 0 by the operators Yt(r) = S^r)CltS{r) taken at r = —to; for all t > 0 
these operators have the shifted form 

Yt(r) = R(r) (g) 1 * -(- i / (g) q* = Rt(r) -f Qt, (2.5) 

K 

where R(r) = ?7l(r)Rt7(r) and Ij is the operator 1 for t > 0 and 0 for t ^ 0 in 
L^(A). It follows from the commutativity condition 

[Y“,Y,^] = 0, Vs,t; a,/3=l,...,d (2.6) 

that the observables {Y(} are selfnondemolition in the sense of their joint mea¬ 
surability, and they are nondemolition with respect to an arbitrary particle 
operator Xt{r) = Ui{ty{X g) l)Ui{t) at to = —f in the Heisenberg picture in 
the sense of their predictability mini 

[Y„Y“] = 0, Vs^t; a=l,...,d; 

indeed, Yg“] = 0 and Yg = Yt for s, t > 0 since 

Ul{s)CltUi{s) = S\-to)QtS{-to ), Vs ^ t. 

Let us fix a state vector /o G A^(A), ||/o|| = 1, given by a localized wave 
function /o(y) on A at y = 0. Let |y) denote the (generalized) eigenfunction 
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in the spectral representation q = /y|y)(y|dA. In the discrete case one can 
take the sharply localized /o = given by the function /o(y) = e(y)/e'^/^, 

where e(y) = 1 if y = 0 and e(y) = 0 if y 7 ^ 0. This yields the localizing 
transformations Ft{y) = (y|5't/o in the form 


Ft{y) = /o(y/- Rt) 


fo{yI - R), t>0 
foi.y)l, t^O. 


The reduced transformations ry 1 —> ipit^y), y G A, defined on the particle 
space Ji by the formula 


i^it,y) = U{t-to)Ft{-to,y)r] , y e A, 

with U{t) = exp{—i/hF[t} and Tt(r) = W{r)FtU{r) reproduce the unitary 
evolution Ui{t) on 77 0 /o G similarly to Eq. (II.dll and namely, 

■*/'(^,y) = U{t){y\StfoU{-to)T]= (y|C/i(i)(»7«'/o), VryeTf. 

The operators Ft{r) at r = —to with an initial wave function 77 of the particle 
before the scattering (tg ^ 0 ) define the probability measure 

ll^(i,y)? 7 fdA = (77,nt(A)77), A C A, 

JA 

for the statistics of the nondemolition measurement of kR via the observation 
of the pointer position y G A after the scattering. It is given by a positive 
operator-valued measure nt(A) = |/o(y — «:Rt(—io))PdA, which is normal¬ 

ized, lit (A) = d, since the measure dA on eA = or A = is translation- 
invariant. 

By renormalizing the operators F(q) as F(y) = e'^/^F(y) with step e = k, 
one obtains the orthogonal projections 

£’(y) = e(/cR - yJ) = / |x)(x|dx. 

Jx:[x/Kj=y/K 


in the case kR = f A(x)|x)(x|dx when A(x) = k[x/kJ . Thus, lit (A), A C A, is 
the spectral measure y) of the quantized position kR(— tg) of the 

particle for a t > 0 given by the eigenorthoprojections E(r, y) = l7^r)E(y)l7(r) 
of the operators R(c), corresponding to the rescaled pointer integer values y/K. 
Thus, the projection reduction postulate has been deduced >from the Hamil¬ 
tonian interaction mi and the nondemolition measurement for the sharply 
localized initial state /g. But there is no continuous limit as k ^ 0 of such 
sharp reduction with nontrivial e 7 ^ 0, since ||e|p = / e(y)dA = ^ 0 and the 

sharp function e(y) disappears as the element of the meter state space T^(A). 

In the continuous case, one can take an unsharp /g G A^ (A) and renormalize 
the operators (??) as Gt(y) = /g(y/ - KRt)//g(y) if /g(y) 7 ^ 0, as in the 
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Gaussian case (II3J. They define Gt(y) as the identity operator for t ^ 0 or 
r > 0, whereas for t > 0 and r ^ 0 the operator 

G(y) = (y|^/o = (y|G, G = f^^Sfo, (2.7) 

say, of the Gaussian form 

G(y) = exp{7rKR(y/ - i kR)} 

given by the generalized eigenfunctions |y) = |y)//o(y) for the spectral rep¬ 
resentation q = f y|y)(y|d/ig with respect to the initial probability measure 
d/j,Q = |/o(y)pdA with density |/o(y)P = exp{—7ry^} in the case 111.511 . The 
corresponding propagators 

T'(t,y) = U(t-to)Gt(-to,y), y e A, to < 0 

define the operator-valued measure nt(A) as 

nt(A)= [ T\t,y)T{t,y)dfio= [ |Gt(-to,y)pd/ro, 

J A J A 

SO that the output probability measure /rt(A) = {r],Ilt{A)r]) is absolutely con¬ 
tinuous with respect to fJ-oiA) = 0 ^J-t{A) = 0. Hence, the reduced 

state vector x(t,y) = T{t,y)r] is normalized to 1 as a stochastic vector process 
x(^) ■ y y) & in the mean square sense with respect to the input 

probability measure /Xg, 

llx(i)llo = J {x{t,y),x{t,y)}df^o = {v,v) = ^- 

This model of nondemolition observation with continuous data y S also 
applies to unsharp measurement of operators R with discrete spectrum. In 
contrast to sharp measurement, unsharp measurement is not sensitive to the 
continuous spectrum limit as k ^ 0 of R = [x/kJ, corresponding to the re¬ 
placement of kR by X G 

Theorem 1 For any initial to ^ 0 the stochastic vector process x(t) = T{t)r] 
satisfies the single-kick equation 

dx(i) + ^Hxit)dt = dlt[G - /]x(t), x(fo) = V, 

generated hy the random differential dlt[G—/](y) = (G(y) —/)dlt ony G A with 
respect to the initial probability measure yiQ. This simplest reduction equation 
is written in terms of the forward differentials dx(t, y) = x(^ + dt, y) — x(^; y); 
that is, is understood in the sense of ltd. 
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Proof. Indeed, by representing G* as Gt(y) = (G(y) —= J + lt[G — 
/](y), for x{t) = U{t — to)Gt{—to)r] one obtains the integral equation 

X(<) = U{t- to)r] + U{t)lt [G - I]U{-to)r) = 



But this equation is equivalent to the differential equation (??), a fact that 
can be proved by straightforward differentiation taking into account the Ito 
multiplication table 


(dt)^ = 0 , dMlj = 0 = dljdt, (dl*)^ = dl^. 


Similarly, one can obtain the simplest nonlinear stochastic equation for the 
normalized reduced state vector Xy{t) = x(^iy)/llx(^:y)ll- 


<^Xy{t) + '^/fiHXyi.t)dt = {Gy{t) - I)Xy{t)dlt , Xy(0) = r], 


where Gy{t) = G(y)/||G(y)Xj^(t)||, rj € H. This equation is an equivalent 
differential form of the nonlinear integral stochastic equation 



= U{t - to)v + U{t){Gy - /)ltx,(0) = C/(t)Gt(y)x,(0)/||Gt(y)x,(0)||. 

This yields Xy{t) = T{t,y)'q/\\T{t,y)-q\\ for an to ^ 0 because of ||G4(y)7?|| = 
\\T{t,y)T]\\. 


Note that the random state vector Xy (t) is obtained by conditioning with re¬ 
spect to the output (rather than input) probability measure d/i = ||x(i; y)lPd/io. 

3 Spontaneous Localization of a Single Particle 

Let us consider a spontaneous process of scattering interactions HU of a quan¬ 
tum particle at random time instants tn > 0, ti < t2 < ■ ■ ■, with a renewable 
meter in an apparatus of the cloud chamber type with bubbles serving as the 
meter. We consider the increasing sequences as countable subsets 

r C R+ such that Tt = r n [0, t) is finite for any t ^ 0 in accordance with the 
finiteness of the number of scattered bubbles on the finite observation interval 
[0,t). The set of all such infinite r will be denoted by Too, T is the inductive 
limit UTt as t > oo of Tj = {rt : r € Too}, which is equal to the disjoint union 
r* = J2'^=o^tin) of n-simplices rt(n) = {G < • • • < c 

The measurement apparatus is assumed to be a quantum system of infinitely 
many bubbles each of which is identical to the single meter described in the 
previous section. The coordinates a = 1,..., d, of a bubble labeled by the 
scattering number n G N show the position y„ G A of the pointer at time tn G r. 
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The interaction Hamiltonian of the particle corresponding to these scatter¬ 
ings is given by the series 


i?(<, r) = iJo + kR 0 ^ (5(t — t„)P(n) (3.1) 

n—1 

having at most two nonzero terms when t G r. Here Hq = iJ 0 1 is the Hamil¬ 
tonian describing the time evolution on the intervals between the scatterings 
t £ T and P(n) is the momentum of the nth scattered bubble, given as the 
vector-operator P(n) = {Pi{n),..., Pd{n)), where Pa{n) = —ihA/Aq’^ in the 
case A = and R 0 P(n) = ^ Pa{n)R°‘. 

The generalized Schrodinger equation corresponding to the Hamiltonian 
ED can be written for fixed r G Too by analogy with the single-kick case 

d'(/'(i) + ^ RoV'(i) dt = (5'(nt) - J)'0(t)dnt, V'(0 ,t)=V’o- (3-2) 

Here S'(n) = exp{—(i/?i)nR (g) P(n)} and nt{t) = \Tt\ is the numerical process 
that gives the cardinality |Ti| = the localized subset < t}, 

so that Ant{t) is equal to 1 for t G r, and zero otherwise. 

Proposition 2 The solution to the equation ‘A) is uniquely determined for 
every r G Too by the initial state ipQ of the system. Namely, -ipftjT) = U{t,T)ifQ, 
where U{t,T) = Uo{t)Vf (t), {t) is the chronological product = 

... S{ti), and 


Vt{T) = sUh)SUt2)--- 



(3.3) 


Here St{tn) = ffg (tn)<S't(n)17o(tn) for tn < t, where S(n) = exp{—(i//i) kR 0 
P(n)}, and Stftn) = I for tn ^ t, so that the infinite product id.di) contains only 
a finite number Ut = X^rGr ^t-r of factors different from the identity operator 

I. 


Proof. Recall that the differential equation (E3 is equivalent to the integral 
equation given by the recurrence relation 


i;{t, t) = U, + - I^r, r)] (3.4) 

\ rG T ^ 


for every r G Too. Hence, 'f{t,T) = C/o(t)V)t(r)^Q, where Uo{t) = e '^ot/h 
Vt(T) is a solution to the operator equation 

r<t 

Vtir) =I + J2 ^r{t){Sl{r, t)-I), Vo{t) = / , 

r^T 
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where S^t,T) = Uo{tyS{nt{T)yUo{t). But this equation has a unique solution 
(ira . which can be written as the binomial sum 

[Ltih) + I][Lt{t2) +/]•••= ^ L{ri,T) .. .L(r„,r) 

pCrt 


in terms of p = {ri,.. . ,r„}, ri < ■ ■ ■ < Vn, n ^ m, Lt{r) = Sl{r) - / (= 0 
if r > t) and L{r,T) = S''^{r,T) — I. Indeed, this sum contains I as the null 
product corresponding to r = 0 and the sum of the other terms is equal to 

r<t p<r 

VtiT) - I = ^ L{ri, r) . .. L{rm, T)L{r, r) 

rGT pCTr 

r<t r<t 

= Vr{T)L{r, r) = ^ K(r)(S'^(r, r) - J), 


where m ^ nj — 1. ■ 

Note that the differential equation (E3 depending on r G Too via n* = nt(r) 
is not stochastic as long as we have not fixed a probability distribution for the 
instants r = . ■.) of the spontaneous interactions. To obtain a continuous 

(at least in the mean) dynamics for such an instantaneous process, one can 
assume that the probability distribution of the random number process ntij) 
is given by the Poisson law 7ro(dT) on Poo presented as the projective limit as 
t ^ oo of the probability measures 

’’■o(dTt) = dr^, v>Q. (3.5) 

Here rt = r is a finite time-ordered sequence r(n) = (ti,...,t„) S P* with 
n = rit, dr* = rifcli dtfc is the measure on Ft given by the sum of product 
measures dti,..., dt„ = dT(n) on the simplices rt(n), dr(0) = 1 on rt(0) = {0} 
such that 


: = 


n—0 


dti... dtn = e’' 




Note that any other numerical process can be described by a positive density 
function /(r) with respect to the Poissonian measure, that is, has the form 
/(r)7r(dr). 

Let us fix an initial state Pq = /“ of the apparatus as the infinite product 
fo° = of the identical state vectors fk = fo of the bubbles given by a 

normalized element fo G L^(A). We suppose, as in §1, that /o(y) 7 ^ 0 for almost 
all y 7 ^ A such that the space £ = {p : A ^ C : ||(/3||o < oo} is isomorphic to 
L^(A) with respect to the multiplication p{y) = /(y)//o(y) by and the 
scalar product ||(^||q = / \p{y)\^dpo = II/IP- This defines the solutions 
of the stochastic equation lira with the initial data ipo = given by state 

vectors rj G H in the particle space Ti. as the state vectors in the product space 
Tdao —'H® £oo, where £oo = lim„^oo L^(A"’) is the Hilbert space generated by 
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the infinite-product functions ^p{y) — O^i fkijk) with fk = fo for almost all k. 
One can describe £oo as the space of all functions (p : A°° —> C square integrable 
with respect to the product measure ^g°(dy) = ^Q(dyi)-/rg(dy 2 ) • • • on the space 
A~ = AxAx... of sequences y = (yi,y 2 ,...): Ml = / < oo). 

The generalized product vectors \y) = |yi) 0 |y 2 ) 0 ... of this space are defined 
by the ^-functions (q|y) normalized with respect to Pq : f |y)(y|yg°(dy) = 1. 

Consider the sequence q = (qi, q 2 , ■ ■ ■) of coordinates (/“, a = 1,... of 
the scattered bubbles at the time instants {ti,t 2 , ■ ■ ■ }■ The commuting vector- 
operators q„, n S N, described in £oo by the multiplications qn|y) = \y)yn, 
are assumed to be measured at the random time instants tn, n G N. The point 
trajectories of such measurements are given by the sequences y = (yi, 2 / 2 , • ■ •) 
of pairs y„ = (tn,yn) with ti < t 2 , - ■ ■ and y„ G A, identified with countable 
subsets V = {yi,y 2 , ■ • ■} C M+ x A. As elements v = {T,y) of the Cartesian 
product Too = Too x A°°, they have the probability distribution Po(du) = 
7ro(dr)/rQ“(dy), where A“ is the space of all sequences y = (yi, y 2 , ■ ■ ■ ),'yn S A, 
equipped with the probability product-measure /ig°(dy). 

The measurement data of the observable process up to a given time instant 
t > 0 is described by a finite sequence Vt = (yi,... ,y„) with n = nt{v) given 
by the numerical process ntir) for the component r of v. 

Let us introduce the counting distribution nt(A) = \vt C (M+ x A)| as the 
number nt{A., v) of scatterings in the time-space region [0, t) x A and define the 
counting integral L(r,y)dnr(dy) over y G M+ x A as the series 

n[L](u) = ^L(y), VuGToo- (3.6) 

y^v 

Having fixed an integer-valued distribution nt{A) G {0,1,...} as a function of 
t ^ 0 and of measurable sets A C A, one can obtain the corresponding trajectory 
u as a sequence of the counts of the jumps of nt{A) in the time-space R+ x A. 

Given an initial state vector in Tiao of the form tpQ = y (g) Pq with fixed 
ifQ ~ /“, one can define a nonunitary stochastic evolution rj i-^- T{t,v)rj by 
setting 

T{t,v) = {y\U{t,T)ipf^, v = {T,y), 

which reproduces the unitary evolution U{t,T) = Uo{t)Vf (t) defined by 13.211 . 
This can also be written as T{t,v) = U{t)F^{v), since Uo{t) = U(t) ® 1 com¬ 
mutes with the (generalized) eigen-bras (y| = (yi,y 2 , ■. ■ | of the bubble coor¬ 
dinates (qi,q 2 ,...) : (y|C/o(t) = U{t){y\. The reduction transformations Ft{v), 
V G Too, are given by the chronological products 

Ft{v) = Gl{y^)Gl{y2) ••• = [] Gl(y) (3.7) 

yGv 

of Gt{tn,yn) = U^tn)G{yn)U{tn) for tn < t, where G(y) = {y\Sfo, owing 
to the product form lES of the unitary transformations Vt{T), r G T, and 
{y\ = Ofci(yfe|> V’oiy) - O^i /o(yfc), {yM = 1 for y e A°“. The stochastic 
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operator (IT7I) defined by the single-point reductions 


Gt(r,y) = {y\St{r)fo 


G{r,y), 

J, r > t 


(3.8) 


is normalized with respect to the initial probability Po(dr;), nT^[/](t) = I, 
where 

n^[X](t)= [ T\t,v)XT{t,v)PQ{dv) = [ Ft{v)XFHv)Po{dv), (3.9) 
J A J A 

is a continual operational-valued measure 0-0 defined on measurable sets 
A C Too of the point trajectories Vt = {(r, y) G v\r < t} given by the operations 
<i>t(r;) : X Ft{v)XF^ (v) for particle operators X : H ^ H. 

The positive operator-valued measure Hi (A) = nyi[/](t) gives the statistics 


Pt{A)= j ||(y|f7(t,r)7/>o||Vo (dy)7rt(dT) 

JA 

for the continual observation with respect to an arbitrary initial wave-function 
= r] ® rj ^Ti, VA the form 

Pt(A) = (7?,nt(A)?7). 

The output probability measure P(Gl), A C Too, is defined by the marginales 
Pt(A), A C Tt, as t —> oo. 

Theorem 2 The reduced wave function x(t,t;) = T{t,v)r] is normalized 


llxWf = J llx(i>'(^)f Po(du) = 1 


as a stochastic vector process xif) : Too ^ Tl with respect to the initial proba¬ 
bility Pq. It satisfies the stochastic wave equation 

dx(0 + ^^X(0dt = dnt[G'-/]x(t), x(0) = ry, (3.10) 

expressed in terms of the random differential dntjG — I]{v) = (G(y„j(u)) — 
l)dnt(t;), nt{v) = nfir) for the point distribution nt{L\ = /^L(y)nt(dy) = 
n[Lt] over y G A defined as hH.bt with Lt{r,y) = \t{r)L{y). 

Proof. To prove Eq. (IXTUIl . discovered for the first time in m , we rewrite 
it in the following integral form 

X(t) = / X " /)x(r)dn,(dy)^ , 

given for each v G Too by the finite sum 

r<t 

n[ltU\r){G - I)x]{v) = ^ {r){G{y) - I)x{r) ■ 

(r,y)Gv 
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We express the solution to this equation in the form x(t, f) = U{t)F}{v)'q via 
the solution to the recursion equation 

r<t 

Ft{v)=I+ Fr{v){G{r,y) - I), Fo(u) =/, 

{r,y)Gv 

with G{t,y) = W{t)G{y)U{t), as was done for the unitary case. ■ 

Let us also write on F the nonlinear equation 

^ = (G„(t) - I)x^{t)dnt{v), x,;(0) = V, 

with Gv{t) = G{yrn{v))/\\G{ynt{v))Xvii)\\- Ks solutions define the normalized 
reduction x^(t) = x(t, u)/||x(L 1^)11 for continual counting measurements as a 
stochastic vector process Xvi^) G ^ respect to the output probability 
measure P of the point process t vt- This can easily be obtained, as in d, 
by applying the ltd multiplication table 

(dt)^ = 0 , dt dnt = 0 = dritdf , (dnt)^ = drit. 


4 Mixing Reduction for Many Particles 

We now consider M identical particles interacting independently with the bub¬ 
bles in accordance with the scattering term in the Hamiltonian (Id. HI . The 
spontaneous process of scatterings is described by the time-ordered sequences 
of pairs (fc„, t„), ti < t 2 < ..., where G M} is the number of the 

particle labeled by the scattering number n G N at time instant > 0. We have 
excluded the possibility of two or more scatterings of the bubbles at the same 
instant of time, as was done for a single particle in §2. The sequence (ki,ti), 
(^ 2 ,^ 2 ),... of the scatterings can be represented by the occupational subsets 
Tk = {tn G T : = fc} of the time set r = {ti,t 2 , ■ ■ ■}, which are disjoint, 

Tfc n T; = 0 if fc 7 ^ Z, since the scatterings for different particles are independent. 
We shall consider the M-tuples t, = {ti, ... ,tm) of these countable subsets 
Tk C R+ as elements r, G T^ of the Cartesian M-product Too, given by the 
partition r = Ut^ : = Ur^, (b t; = 0 if fc 7 ^ / of a r G Too. 

The interaction Hamiltonian of M particles with independent scatterings 
labled by a sequence t, = {(fci, fi), (^ 2 , ^ 2 ); • ■ • } reads 

00 

H{t,T,) = -I- R(A;„)) 0 5{t - t„)P(n). (4.1) 

n—1 

Here = Fl^ 0 1 is the Hamiltonian of the particles describing the time 
evolution on the intervals between the scatterings with the bubbles: 

M MM 

W{k,l), 

k—1 k—1l>k 
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where H{k) = is the Hamiltonian of the /cth particle 

and W{k,l) is the interaction potential in of the fcth and Hh particle, 

I ^ k <l ^ M. 

Let denote the M-particle Hilbert space, which is an invariant subspace 
■j^M Q pj: symmetric (bosons) or antisymmetric (fermions) M-tensors 

r]^ € generated by the product-vectors G with t]^. € Ti.. The 

correspondent Ito-Schrodinger equation for the stochastic state vector (i) ■ 
T, 1-^ T.) with values r.) G Ti.^ in the product space Ti^ = Ti.^ 0 £oo 

of the M-particle space by £oo = lim„^oo L^(A”) reads 

^ = {S{kt, nt) - {t)dnt ■ (4.2) 

Here S{k,n) = exp{-i kR(A:) 0 P(n)}, R{k) = 0 R 0 

kt{T,) ='^klr^{t), where l^(t) = i ^ 

fc=l [U, tfT 

is the random number kt : T^ ^ {1,...,M}, labeling a particle by k at any 
instant t G of its collision with a bubble labeled by ntij) = J2^=i‘'^k,t = 
\t n (0, t)\, where nk,t = hk n [0, t)|, t = Ur^. 

Proposition 3 The solutions ^{t,T,) = U{t,T,)'tjj^, G of Eq. y.g| ) 

can be written as U{t,T,) = {t)V^ {t,) in terms of the finite chronological 

product 

Vt{T,) = S\{k^M)sKk2M)---, T. Gr“, (4.3) 

where St{kfin) = U^\tn)S{k,n)U^{tn) and U^(t) = ex.p{-j^H^(t)}. 

The proof is exactly the same as for the case of a single particle (M =1). 

Let uj = {wi,W 2 , ■ ■ ■) denote a chronologically ordered sequence of triples 
Wn = ikn,tn,yn) and fl the space of such sequences with {ti,t 2 , ■ ■ ■} & Loo. 
Every sequence w G H can be represented as a pair uj = (T,,y), where r, = 
(ri,..., tm) is a partition of the corresponding sequence r = {ti, ^ 2 , ■ • ■ } and 
y = (yi) y 2 , • ■ • )> so that can be identified with the product T^ x A°“. The 
space H is equipped with the probability measure Po(da;) = 7ro(dT,)/ig°(dy), 
where 7ro(dri,..., drM) = O^i 7''o(dTfc) is the product of the identical Poisson 
measures 113.5ll . in accordance with the independence of spontaneous interactions 
of each particle with the bubbles. 

Given an initial state vector G ipg, where (pg ~ /^, one can easily 

prove that the nonunitary stochastic evolution 

T{t,uj) = {y\U{t,T,)(pQ, uj = {T,,y) 
is also a finite chronological product 

T{t,u;) = U^it)FHu;), Ft{u;) := GI{w,)GI{w2) ... . 
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Here Gt{k,tn,y) ■= I for tn > t and 


Gt{k,tn,y) := U^\tn)G{k,y)U^(t^), ^ f, 

where U^{t) = exp{— H^t}, is defined by the reduced scattering operator 

G{k, y) = 0 G(y) 0 , G(y) = {y\Sfo (4.4) 

for fo G 47^(A), ||/o|| = 1, applied to the fcth particle only in Ti.^. The stochastic 
operator T{t) defines the solutions x^(t,a;) = T{t,uj)r]^ to the Ito differential 
equation 


M 

dx^it) + ^H^x^it) drikAGik) - I]x^{t), x^ A) = 

k=l 

for the stochastic vector states H ^ of the M-particle system, 

corresponding to an initial rj^ G Ti.^. The right-hand side of this equation is 
written as a point integral nt[L] = 'Yk=i f ^(k,y) drit = n[Lt] with respect to 
the stochastic distribution n[L](a;) = L{w) forL(fc,r, y) = lt(r)L(fc,y) = 

Lt{w). The vector x^(t, w) G TL^ as well as x(t,u) in 13.1 OH . is no longer 
normalized (||x“(t,a;)|j A 1) for an initial state-vector rj^ G , Wv^W = !> 
but it is normalized with respect to the probability measure Pq on H in the 
mean square sense. But, in contrast to xi^^A^ vector x*^(i,w) is not yet 
the reduced description of the M-particle system under the observation of the 
scattering process Vt = {{r,y) £ v : r ^ t}, given by the registration of the 
pointer positions yn G A at random time instants tn- 

The reduced dynamics corresponding to the observation is described by a 
stochastic operational process X ^ 0[A'](t), 

0[X](t,r;) = <^t[U^At)XU^ {t)]iv) (4.5) 

for the M-particle operators X : Ti^ —> Ti^ given by the conditional expecta¬ 
tion 

$t[X](T,y) =Y Ft{<y.,y)XFAcr„y) (4.6) 

Uo-fc=rt 

where the sum is taken over all partitions cr, = (ui,..., (Tm) of a finite subset 
Tt = tC] [0, t). This averaging is due to the impossibility to detect the individu¬ 
ality of the identical particles producing indistinguishable effects on the bubbles 
by measuring the scatterings of the bubbles. 

To prove Eq. (lOll . we need to compare the correlations of Ft{Lo)XFl{uj) 
and of an arbitrary functional g{vt) of the observable point process Vt with the 
correlations of 14.611 and of g(yt)- But by applying the well-known formula pT] 
one can easily find 


M 


c(cri,...,f7M) dcT/e = 




4 S = 

^ t (Tfe:U(Tfc=(T 


:(tTi,...,crM)dCT 
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for the multiple point integration that these correlations with respect to the 
probability measure Pq on fl given by the Poisson law (1^ simply coincide: 



M 


tM 


t k—1 


f (ff(CT), X! ^(CTi,...,crM))oe 

/T>.'\ l/Tl_=/T 



Tk-'JTk—O' 


Here (•, •)o is the abbreviation for the inner product in £ of the test func¬ 
tion V I—> gir^y) with fixed r G P* and the operator function Xt{Tt,y) = 
Ft{T,,y)XF^{T,,y) with fixed r, S P^ and r = U^^Tfe. The probability 


measure 


M 



UCTfc=rt k=l 


on Pt has the intensity Mi/. It is induced by the measure 7ro(dr,) on P^ with 
respect to the particle identification map r, € P^ t = Uk=iTk, defining 
the observable data v = {T,y) by the stochastic map co = {T,,y) (r = 
y) G Too on w G H. 

By the coincidence of the correlations proved above, the stochastic operator 
(14.511 is indeed the conditional expectation of the stochastic operators Xt{uj) 
with respect to the observable process t = Vf 

In contrast to the pure operations X i—> Xt{uj), the reduction operation 
X I—> $t[X](z;) preserves the symmetry of the M-particle operators X with 
respect to particle permutation. It is the least mixing operation which preserves 
the indistinguishability of the particles with respect to the observations of the 
bubble scatterings, corresponding to the complete nondemolition measurement 
of the particles. 

Indeed, the reduction operation (14.(ill can be simply written as the finite 
iteration 



= ^'t[... 4't[X](yi).. .](2/n) = ...,yn) 

with n = \Tt\ single mixing reductions 



(4.8) 
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Given as the arithmetric mean value of the permutations for the pure operations 
X I—> G\(k,y)XGt{k,y), corresponding to the identical operators (14.411 . the 
reductions X i—> 4('i[X](y) are permutationally symmetric and are not mixing 
only if the pure operations do not break this symmetry. 

The mixing property of the reduced stochastic dynamics t p[X]{t,v) de¬ 
rived above for the corresponding statistical states 

p^[X]{t,v) = = Tr{Xp“(t,u)} (4.9) 

gives an increase in the entropy 

a’^{t,v) = -Tr{p^(t,t;)lnp“(t,u)} 

for an ensemble of identical particles even under the condition of complete 
nondemolition observation. According to (lOll . the reduced density operators 
p^{t,v) for the system of M identical particles gives the probability density 

= TT{p^{t,vt)} = P^{dvt)/P^{dvt) 

of the output process Vt- Here G is the probability measure on 

Too = Poo X A°° defined by the Poisson measure SII). This means that the a 
posteriori density operator p^ {t) is defined as a stochastic positive trace class 
operator normalized in the mean sense 

\\p^{t)h= I Tr{p^{t,v)}P^{dv) = l. 


Theorem 3 The density p^ (t) satisfies the stochastic operator equation 

, M . 

^ [H, P^{t)\dt = dnJ — ^ G{k)p^{t)G\k) - p^{t) j, (4.10) 


fc=i 


which has a unique solution for every initial condition p^{0,v) = p^ given by 
the density operator p^ for the M-particle states p^[X] = Tr{Ap^}. 

Proof. Let us prove Eq. lUTTTll written in the equivalent integral form 

p^[X]it) = p^[X{t)]{t)+ f f p^mX{t-r)]{y)-X{t-r)]{r)dn,{dy) 

Jo Ja 

for an p^[X] = {rj^, Xrj^), G , where X{f) = (t)^XU^(t), 

M 

U^{t) = vl/[X](y) = _^Gt(fc,y)AG(fc,y). 

^ k^l 

Taking into account the fact that the stochastic integral in this equation 
is simply a finite sum for every v G Too and t, one can write it as a recursive 
operator equation 

r<t 

$*[A](r;) = A+ ^ [vf-[A] (r, y) - A](r;), 

(r,y)ev 
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for a stochastic operation defining the solutions to Eq. (I4.1()ll . as in 114. 9II . in 
terms of the composition (Oil . 

But such recursive equation has a unique solution = ^t(yi) o '^t{y 2 ) ° 

..., defined for a t> = (r, y) G Too as the chronological composition of the maps 

if r ^ t and y) = X if r > t. This solution 

can be found by the iterations 


Mv)-I = 


r<t 

E 

{r,y)Gv 


<^r{v)oA{r,y) = 


(r,y)evt 


A(r,y) + 


(s,y)Gv, 


d>s(u)oA(s,y) 


where A(j/) = 'f'(y) — I, I is the identical map X X, and by the binomial 
formula o ^'(( 2 / 2 ) o • • • = o ’ ’ ’ o A(z„) in terms of ct = 

{zi,.. .,Zn}, z = (r,y), ri < ■ ■ ■ < r„, n < n*. ■ 

Let us also write the nonlinear stochastic equation 


dpf (<) + ^ [H, dt = p^(t) o - I) dnt(v), 

where = ^'(y„^(„))/Tr{E(y„^(„))p}, 

M 1 

po^(y) = M 

fc=i fc=i 

for the normalized density operator [t) = p^{t,v)/p^{t,v). This describes 
the conditional expectations p^[X]{t) = Tr{Xp^(t)} of the M-particle opera¬ 
tors with respect to the output probability measure P^fdu) = p'^ft, t;)P/}^(du), 
where p^(t, u) = Tr{p“(t, u)}. 


5 Macroscopic and Central Limits of the Model 

We now consider the mean field approximation of the measurement apparatus 
fixing its total effect vk = —7 given by the mean number v of scattered bubbles 
per second and an interaction constant k coupling each bubble to a particle 
in the Hamiltonian eu. We look for the limits of the unitary and reduced 
evolutions (ESI and lirnni as v 00 and k —> 0 such that 7 is a real constant. 
To perform these limit passages, we need the expansions 

S'(n) = 7 (g) 1 - i^ R (g) P(n) “ Q ^ ® P(^))^ + • ■ • 

of the scattering operator S{n) = exp{—KRgP(n)} and the reduced operator 
G'(y) = fo{yI — 'iR-)//o(y) with respect to the coupling constant k. The first 
term of the expansion for S{n) disappears in the right-hand side in Eq. lid. 211 . 
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whereas the second and third terms give rise to the differentials of the operator¬ 
valued stochastic integrals 


/o 

The corresponding terms 

'■ m) 


= [ 'Pinr)dnr, nJP^] = f P(n^)^dn^. 
Jo Jo 


Jo 


0 JA 


/o(y) 


dnr(dy), nt[ ^ ] = 


0 Ja 


/o (y) 

/o(y) 


dnr(dy) 


on the right-hand side in Eq. can also be written as the integrals fit [L] = 

Jg L(rir) dur with values in operator functions of t € Too 


nt(T) 

fit[L]{T) = Y, L{n), Lin) = ® L 0 , (5.2) 

n—1 


where L is one of the multiplication operators 

i' = /^(q)//o(q), L" = /^'(q)//o(q). 

Here /^(q) denotes the vector of the derivatives 9a/o(q) = dfo{<^/dq°‘ and 
/"(q) denotes the matrix with elements dadfif{q). The operator nt[L\{T) cor¬ 
responding to a function L = Z(q) of the bubble coordinate vector-operator q 
in T^(A) acts in Eoo as the multiplication operator 

nt{T) 

[nt[L]{T)(p]iy) = Y ^iyn)y^iy) = nt[l]iT,v)ip{y). 

n—1 

Hence, the main terms on the right-hand side in Eq. (E3 and (liriT)|l for 
K —*■ 0 are given by the renormalized stochastic integrals 

A(f) = ^ / Liur) dur = \ nt[L] (5.3) 

nt Jq ut 

of the operator-valued stochastic functions L{t,T) = Lintir)) with respect to 
the numerical process ntir) that has the Poisson probability distribution (13.511 
on Too- 

To pass to the large number limit > oo in itCT for an arbitrary oper¬ 
ator L in T^(A), we need to use the quantum stochastic representation |22 
of the integral 15.411 in the Fock space F over L^(R+ x A). The space F 
is defined as the L^(T)-space of all square integrable functions (j) : T ^ C, 
II^IP = /y |(/)(r;)pA(du) < oo of time ordered finite sequences v = {yi,... ,yn), 
y = {t, y) identified with subsets v C R+ x A of finite cardinality |u| = 0,1, 2,.... 
The measure A(dt;) on the union 


T = Y 

n—0 
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of the disjoint subsets T(n) = {u e T : |u| = n} is given as the sum 

CXD 

A(A) = ^A(T(n)nA) 

n—0 

of the product A(du) = measures dy = dt dA on R-|_ x A such that 


= E 

n—O 




i=l 


Let Nt[L] denote the numerical integral as an operator 


nt(T) 

(Nt[L](f>)(T,y) = ^ L(n)ip(T,y) = nt[L](r)(/)(r,y) (5.4) 

n—1 

in the Fock space (j) G L^(T). This integral represents the operator-valued 
stochastic integral itCTi by pointwise multiplication of the functions (/>(v) = 
(l>{t,y) of u S T, by nt[L](T), which is considered as the function of a finite 
sequence r G F because of its independence of T[( = {tn > t}- In order to 
obtain the initial probability measure Po(dt>) = 7r(dr)/rg°(dy) on Too induced 
by an initial Fock vector (pQ G L^(T), we need an isomorphic transformation of 
(15.511 

Nt[L\ = Ni[L] + V^{At[flL] + Al[Lfo]) + (5.5) 

which can be locally performed by a unitary transformation 

N,[L] = UlNt[L]Us, U, = exp{y^;(4[/o] - (A^^])} 

for every t < s. Here A|[/] and At[P] are the creation and annihilation integrals 
of / G L^(A), /I G L'^{A)*, given by the operators 

yevt 

(At[p]cj)){v)= [ f{y)*(l){vUy)dy 

^ ' J[0.t)xA 

in the Fock space L^(T), where v\y means the sequence u G T with deleted 
y = (r, y), r < t, and uUy means the sequence t> G T with an additional element 
y ^ V. The characteristic functional of the stochastic operators nt[L] with 
respect to the initial state-vector ip^ ~ in and the Poisson probability 
measure (13.511 is now given simply by the vacuum expectation 


^^(<^o,e™‘[^Vo)^o(dT) = (<5^,e*^‘W5^), 
where S^{v) = 1 if u = 0; otherwise, S^{v) = 0. 


20 


The corresponding representation l{t) = ^ Nt[L] for 15.411 helps us immedi¬ 
ately obtain the quantum large number limit 

lim \Nt[L] = f^Lfoi 

U^OO l/t 

as the mean value Iq = (/o,T/o) = fo^fo of a single-bubble operator with 
respect to an initial wave packet /o G L^(A). This gives the following: 

Proposition 4 The macroscopic limit 

dt = 4 7(R- ® Poi)V'(^) dt 

n n 

of the generalized Schrodinger equation iTO) turns out to be nonsingular with 
an additional potential —ypoR corresponding to the mean momentum po = 
(/o,P/o) of a bubble in the initial state /o- It coincides with the large number 
limit 1 / —> oo o/ the reduction equation under k = —^fv. 

Proof. Indeed, the mean field dynamics preserves the product structure 
ip(t,v) = r]{t)(pQ{v) of the initial product-vector = g ® on the Fock 

component (pQ G !F because it is identical by virtue of Hq = H (g) 1. But 
unexpectedly (compare with [IHIl the large number limit 

dx(t) + ^ di = -^ 7PoRx(t) dt (5.6) 

of the stochastic nonunitary equation (ITTUIl corresponds to the same unitary 
dynamics r]{t) = U{t)r] = x{t) of the particle state vector if x(0) = p, since 

^ nt ^ (/o, L'fo) = J /o(y)*/o(y) dy = ^ Po . 

■ 

Note that the macroscopic limits for the M-particle system also give essen¬ 
tially the same type of unitary evolutions in the large space 0 IF and in 
the reduced space TL^. To get this correspondence, one has only to replace the 
measure (lint on Too for Eq. (lO l by the product measure tt® on T^ and for 
Eq. 14.1011 by the induced measure 63 on Too, taking Mv instead of v in 15.511 . 
This means that the mixing property of the reduced equation (OUI) vanishes in 
the mean field approximation for the bubble system. 

Let us now pay attention to the fluctuations with respect to the obtained 
large number limits. Such fluctuations might appear for —k = —> 0 in the 

large time scale t ^ 1/k. We can get these fluctuations without rescaling the 
time t if we assume that po = 0 and —k = y/yT?, so that we have to take into 
account also the K^-terms in (EHi- 

It follows from the Fock space representation (15.611 that the quantum central 
limit 

lim Nt[L] = At[f^,L] + AliLfo] 
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exists for any single-bubble operator L with zero mean value {fo,Lfo) = 0. We 
first apply this central limit theorem to the right-hand side in represented 
in (g) as 

"•I® --'I ='■I - Ks) ®'”''1 + - • 

This yields lim ^ Nt[P] = Ut, lim i fVt[P„P^] = (P„/o, Pf 3 fo)t, and 

d-ipit) + Koi^it) dt = y 7(R (g dut)'!/'(t) • (5.7) 

a 

Here = K K = + \ Rcr^R, cr^ = [cr^^], 

Ut = Ht[4P] + AliPfo] = 23?4[P/o] 

is a Fock space representation of the Wiener vector process u : 1Ut with the 
correlations = S^(9a/o, i9^/o), «,/? = defined by the momentum 

operators P due to the quantum stochastic multiplication formula mi 1221 

dM„ dui3 = dAtlf^Pa] dAllP/^fo] = foPaP/sfo dt. 

The central limit equation (Oil for the unitary evolution of the coupled system 
turns out to be a stochastic Schrodinger-Ito equation of diffusive type driven 
by the Wiener process u — {uj : t G K+}. The same conclusion obviously holds 
for the M-particle system driven by M independent Wiener processes Ut(fc) 
identical to Uj, A: = 1,..., M: 


. . M . 

dtp^{t) + K^tp{t) dA = gf ^ R(fc) ® dU((A;) {t) 


k=l 


where + ^ R(fc)cr2R(/c), ® i. 

We now prove that the application of the central limit theorem to the re¬ 
duction equation (lOHl yields an essentially different type of the stochastic 
evolution, originally derived in mi by quantum calculus method. 

Theorem 4 The central limit v = j kY' —> oo of the stochastic wave equation 
TTm has the following nonunitary diffusive type 


dx{t) + Kx{t) dt = 7 Rx(A) dvt • (5.8) 


Here Vj is a complex Wiener vector-process with the Fock-space representation 

vt = Mfo^'] + MWfo] = 5i^l[(wo + wo)/o] -H i3H|[(wo - wo)/o] (5.9) 

given by the complex osmotic velocity wo(y) = 91n/o(y), d = (di,... ,dd) of a 
single bubble, and the operator K is the same as in (SI) . 
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Proof. Indeed, the central limit of the right-hand side in (ITTUI) with K = 
—7/ y/u yields 


limdnt[G' — /] 


7R lim ^ diVt [L'] -h i R lim i diV* [L"]R 


7Rdvt 


1 

2 


Rcr^Rdt, 


R 


where the representation 15.611 is used for 

nt[L'] = tit and nt[L''] = n* 

so that fQdadfsfo = — (9a/o, cI/ 3 / 0 ) are the matrix elements of lim ^ iVt[L"] = 
/qL"/o. The linear stochastic equation 15.911 obtained has a unique solution 
x{t, v) = T{t, v)f] for a given x(0, v) = rj G H which is not normalized ||x(f, i’)|| 7 ^ 
0 for every Wiener trajectory v : t i—>- Vt but is normalized in the mean square 
sense / ||x(t, ■y)|pPo(du) = 1 with respect to the Gaussian probability measure 
Po of u = {v( : t G K+}. The measure Pq is defined by the zero mean values of 
Vt and by the table 

duc dvp = AAt[flL'^]AA\[L'^J q] = dt 

dv*^dv^ = dAt[f^,L2]dAl[L'^fo] = f^L'^L'^fodt 


of commuting multiplications 


dva dv0 = dvjs dho,, dho, dv^ = dv^ dva ■ 

But the reduction noise V( obtained in the Fock space representation does not 
commute with the real Wiener process Ut = u* in 15.811 . 

dh„ dup = dAtlf^] dd4|[P;3/o] = f^L'^Pfsfo dt 

dua dvfj = dAt[f^] dA\ [L'^/o] = f^PaL^fo dt 

if 5a/3(y) = dadi 3 \nfo{y) ^ 0, since 

[dz)„, du^] = 4 [L'„, P/s]fo df = Y (/o, 9a>3fo) dt. 


In the same way one can obtain the continuous reduction equation 


dp^it) + (Kp^it) + p^{t)K^)dt = 



er’^p^^{t)IL{k) dt -|- 7(dW(Rp'^(t) -|- p*^(t)RdW(), 


(5.10) 
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R = -^ Sfeli where the operator K = is the same as in the equation 

for the unitary evolution of the M-particle system coupled with the bubbles. 
The derived stochastic equation for the M-particle density operator p^{t,w) 
normalized in the mean is driven by the complex Wiener process wt = 
having the same multiplication table as the process VM'Vt in (1^ . The diffusive 
type equation 15.1011 as 14.1011 also has the mixing property. This was also 
derived in by operational method. 
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